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We calculate the 2-triplon contribution to the dynamic 
structure factor of the 2-dimensional Shastry-Sutherland 
model, realized in SrCu2(BOs)2, by means of perturbative 
continuous unitary transformations. For realistic parameters 
we find flat bound 2-triplon bands. These bands show large 
weight in the structure factor depending strongly on momen- 
tum. So our findings permit a quantitative understanding of 
high precision inelastic neutron scattering experiments. 

PACS numbers: 75.40.Gb, 75.50.Ee, 75.10.Jm 



Quantum antiferromagnets which do not have a long 
range ordered ground state, so-called spin liquids, con- 
tinue to attract considerable interest. While there 
are many 1-dimensional examples there are only a few 
2-dimensional systems. Of strong recent interest is 
SrCu2 (1303)2 [1], a realization of the 2-dimensional spin 
1/2 Shastry-Sutherland model [2] 



for the dynamic structure factor which is measured by 
inelastic neutron scattering (INS). 
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where J\ and J2 are the intra- and inter-dimer couplings 
as in Fig. 1. Two spins coupled by J\ are referred to as 
dimers. We focus on Ji, J 2 > 0. 

The state |0) with singlets on all dimers is an exact 
eigen-state of H for all values of J\ and J2 [2] . We found 
|0) to be the ground state (singlet-dimer phase) of H for 
x := J2/J1 below w 0.63 [3] while other results indicate 
an instability at slightly higher values of x, for a review 
see Ref. [4]. 

In many papers it has been shown that the magnetic 
properties of SrCu2 (603)2 can be understood well by 
H in (1) in the singlet-dimer phase [4]. Thus the bo- 
rate constitutes a particularly transparent case of a 2- 
dimensional spin liquid. In view of the extensive spec- 
troscopic data on this system quantitative theoretical re- 
sults for spectral densities are highly desirable. But so 
far only the numerical exact diagonalization (ED) for sys- 
tems of 20 or 24 spins was possible [4] . This approach is 
hampered by the finite size in two ways. First, the ener- 
gies of the excited states display strong finite size effects 
since these states are spatially extended, in particular 
the bound states built from two elementary excitations. 
Second, the ED provides only isolated spikes instead of 
continuous distributions. In this article, we remedy these 
drawbacks by making use of recent conceptual progress 
in the method of perturbative continuous unitary trans- 
formations (CUTs) [5,6]. We provide high order results 
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FIG. 1. Shastry-Sutherland model with spin 1/2 on the 
dots; dimers are solid gray lines. The microscopic angles 
and distances apply to SrCu2(BC>3)2 according to Ref. [1]. 
The primitive vectors a and b span the dimer-lattice r e fi (all 
dimers taken as equal) while a and b span the lattice Tab 
distinghuishing horizontal and vertical dimers. 



We briefly review the derivation of perturbative ef- 
fective operators (Hamiltonian and observable) in gen- 
eral and derive the appropriate INS observable for the 
Shastry-Sutherland model. The corresponding spectral 
density, i.e., the dynamic structure factor, is calculated 
via the T = Green function. We focus on the 2-triplon 
part above the flat, featureless 1-triplon band [7-9]. 

The limit of isolated dimers, i.e., x = 0, serves as start- 
ing point of the perturbative analysis. The basic excita- 
tion is given by promoting one singlet to a triplet. The 
next higher excitation are two triplets and so on. Upon 
switching on the inter-dimer coupling (x > 0) the triplets 
acquire a dispersion and become dressed particles, which 
we call triplons [10]. 

Up to a constant, U in (1) counts the number of 
triplons and we define the particle-number operator Q = 
U + 3N/4 (N: number of dimers). The perturbative 
part V in (1) decomposes into ladder operators V — 
T_i + To + Ti, where the index i denotes the number 
of triplons created (destroyed) by Tj . 

The original Hamiltonian is mapped by a perturbative 
CUT [3, 9, 11] to an effective Hamiltonian (in units of Ji) 



h cB (x) = u + J2 xk c N r w 
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where m — (mi,rri2, ■ ■ ■ ,mk),mi G {0, ±1}. This ef- 
fective Hamiltonian conserves the number of triplons: 
[ifeff i Q] = 0. In each order fc, iJ e ff is a sum of vir- 
tual processes T(m) = T mi ■ ■ ■ T mk weighted by rational 
coefficients C. The sum E' is restricted by the triplon- 
conservation condition mi + • • • + m^ = 0. The effec- 
tive Hamiltonian can be decomposed into irreducible n- 
particle operators H n [5] i? c ff = H o + H i + H 2 + 

The matrix elements of the irreducible H n for the 
infinite system can be computed perturbativcly on fi- 
nite clusters due to the linked cluster theorem. For the 
Shastry-Sutherland model Hq is conveniently set to zero; 
Hi and H2 were determined previously [3, 9] to obtain 
the 1- and 2-triplon energies. 

Applying the same transformation as for H other ob- 
servables are also mapped onto their effective counter- 
parts [5] 

00 fe+i 

O cB (x) =J2 xk Yl 2 Cirrr^Oim^) (3a) 

fc=0 i=l \m\=k 

0{m; i) := T mi ■ ■ ■ T m% _ x OT m% ■••T mh , (3b) 

where O is the initial observable. A useful decomposition 
of g reads 



O c s = Od, n 

n=0 d>—n 



(4) 



where d indicates how many particles are created (d > 
0) or destroyed (d < 0) by Od.n whereas n > is the 
minimum number of particles that must be present for 
Od.n to have a non zero action. For T = measurements 
only the Od,o operators matter. 

The energy and momentum resolved n-particle spectral 
density for the operator O is given by 



where Q^ n > is the retarded rt-particle Green functic 
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Since expectation values do not change under unitary 
transformations the Green function (6) is not altered if 
the the effective operators are substituted for the initial 
ones. Using the decompositions of H and of O as well as 
the conservation of triplons by H e s the individual sectors 
of different triplon numbers can be analyzed separately. 

We focus on the 2-triplon sector and introduce the 2- 
triplon momentum states [3, 12] 



\a, K, d) s = J= ]T e *(K+-QHr+d/2) | F) r + d) f 



(7) 



where S G {0, 1,2} is the total spin of the two triplons 
and |r, r + d) is the state of one triplon on the dimer at 



r and the other one on the dimer at r + d. The primitive 
vectors a and b span the dimer-lattice T e g (Fig. 1). The 
vector K + crQ lies within the (first) Brillouin zone (BZ) 
of the dual lattice T* ff spanned by the vectors a* and b* ; 
as usual Q = (it, n) in units of the inverse dimer-lattice 
constant. The additional quantum number a G {0, 1} 
is chosen such that K lies within the magnetic Brillouin 
zone (MBZ) which is the (first) Brillouin zone of the dual 
lattice r^ B . The exchange parity of the two triplons is 
fixed by \a, K, d) s — (— l) 5 |cr, K, — d) s , hence we restrict 
to d = (di , d 2 ) > [di > or {d 1 = and d 2 > 0)] . 

For fixed total momentum K, Hi (15 th order) is a semi- 
infinite matrix in d and a while H 2 (14 th order) is rep- 
resented by a 84x84 matrix. The matrix elements are 
polynomials in x calculated previously [3, 12]. 

We turn to analyzing the appropriate observable for 
the INS experiment on the SrCu2 (603)2- It reads 
H<1) = EiS^x;)^'* 4 , where S z {xi) is the z- 
component of the spin at the position (dots in Fig. 1). 
The Xj must not be confused with the vectors r G r e ff 
which denote the positions of the dimer centers. 

The momentum transfer q measured in experiment 
is any vector in the dual space whereas the excitations 
of the Shastry-Sutherland model arc labeled best by 
the momenta K G MBZ. The usual backfolding im- 
plies K(q) = q mod (r^ B ) and K(q) + er(q)Q = q 
mod (r* ff ) whence <r(q) = for q G MBZ mod (r* ff ) 
and <r(q) = 1 otherwise. 

We construct an operator Af defined for K and r such, 
that A/"(K;r) and JF(q;x) have the same action on the 
ground state |0). The operator Af will then be used to 
obtain the effective operator. It is a crucial feature par- 
ticular to the Shastry-Sutherland model that the triplon 
vacuum |0) is not changed by the CUT since it is an ex- 
act eigen-state. With a suitable convention for the singlet 
orientation the action of T on |0) is 



S< n >(w,K) = - 7 r- 1 Img(")(o;,K) , (5) ^(q)|0) 
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The sums run over all vertical dimers r v and horizontal 
dimers r h . A state |r) is defined by one triplon with 
S z = on the dimer at r and singlets elsewhere. The 
vectors arc defined in Fig. 1. 

The appropriate local operator W(r) using the dis- 
tances r G r c ff reads Af(r) = Sq(t) — Sf(r), where the 
subscripts and 1 distinguish the two spins on dimer r 
such that Af(r)\0) = |r). The momentum space repre- 
sentation is given by (a = 1 — a) 



A/"(q) = a(q)jV(<7, K) + b{q)Af(a, K) 



<r(q),K(q) 



o(q) = i [sin (q • 8 V ) + sin (q • S h )] /2 
6(q) = i [sin (q • 5 V ) - sin (q • 5 h )] /2 , 



(9a) 

(9b) 

(9c) 
(9d) 
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which ensures 7V(q)|0) = JF(q)|0) for all q. 

Including the microscopic details for SrCu2 (1303)2 
(Fig. 1) and denoting q by q = ha* + kh* fixes the scalar 
products toq-<5„ = 0.717(/i-fc) andq-5 ft = 0.717(/i + fc). 
This completes the derivation of the excitation operator 
A/"(q) for the INS with momentum transfer q. 

The action of the effective local operator A/" e ff(r) from 
Eq. (3) on |0) is implemented on a computer. Although 
Af(r) exclusively produces 1-triplon states when acting 
on |0) A/" c ff(r) leads to states containing an arbitrary num- 
ber of triplons. Focusing here on the 2-triplon channel 
we have to deal with A/2.0 ■ The calculations of the am- 
plitudes of A/2,0 can be performed on finite clusters, see 
Refs. [6,12]. 

By substituting A/2,0 ( r ) for A"(r) in Eq. (9b) we ob- 
tain A/2,0 (q) which excites the same type of 2-triplon 
momentum states \<r, K(q),d) that are used for the ef- 
fective Hamiltonian. The 2-triplon amplitudes A CT ,K,d = 
(a, K, d I A/2,0 (q) 1 0) with K € MBZ define a vector 'in the 
quantum numbers a and d. Each component is calcu- 
lated to 8 th order in x. 

The 2-triplon energy and momentum resolved spec- 
tral density of A" is obtained by evaluating the 2-particle 
Green function (6) via tridiagonalization [13, 14] 



E CT , d K,K,d(x)| s 



a 



ai 



b\ 



(10) 



For fixed K, the continued fraction coefficients <Zj and bi 
are obtained by repeated application of Hi + H2 (ma- 
trix in d and a) on the initial 2-particle momentum state 
|/o) = A/2,0 (q)|0) (vector in d and a). Prior to the eval- 
uation we extrapolate the matrix elements of H\ and 
i?2 and the amplitudes of |/o) by optimized perturbation 
theory (OPT) introduced in Ref. [6]. The values used 
for the OPT parameter «opt are —0.20 for the elements 
of Hi, 0.80 for the elements of H2, and —0.25 for the 
amplitudes of M [12]. 

In Fig. 2, the spectral densities of A/2,0 are plotted for 
the two momenta K = (0,0) and (0,7r), which translate 
to q = (h, k) as indicated. Results are shown for two sets 
of parameters. The set x = 0.635, Ji = 7.33 meV was 
proposed from the analysis of the magnetic susceptibil- 
ity x(T) [15]. We proposed the set x = 0.603, Ji = 6.16 
meV previously [3] based on the analysis of excitations 
energies at T — 0. Clearly, the differences between the 
two sets matter. Both the positions and the weights of 
the curves differ from one set to the other. Comparing to 
high resolution INS data [16] we come to the conclusion 
that the parameters x = 0.603, Ji = 6.16 meV fit signifi- 
cantly better, both concerning the positions and and the 
weights of the peaks. So we favor this set of parameters. 
It is objected that x(T) is not well described [4]. But 
x(T) is also strongly influenced by the presence of inter- 
layer coupling [3,15] which is not known. So a definite 
conclusion on the basis of x(T) alone is very difficult. On 



the contrary, ED data for the specific heat [15] indicates 
that lower values of x fit better to experiment than larger 
ones, cf. the data for 16 spins. 

A possible weakness in our analysis are the necessary 
extrapolations. In Ref. [3], we did not use OPT but Dlog 
Pade approximants which allow for power-law singular- 
ities. This led in a very robust way to the instability 
at x « 0.63. OPT does not allow for power-law sin- 
gularities. Hence it leads to a smoother dependence on 
x. No instability occurs below 0.7 so that the precise 
position of the instability is still an open issue. We em- 
phasize that in spite of the smoother OPT extrapolation 
the parameters x = 0.603, Ji = 6.16 meV still yield a 
better agreement with experiment than the parameters 
x = 0.635, Ji = 7.33 meV. 
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q=(2.0,0.0), x=0.603, J,=6.160meV 

q=(2.5,0.5), x=0.603, J,=6.160meV 

q=(2.0,0.0), x=0.635, J,=7.325meV 

q=(2.5,0.5), x=0.635, J,=7.325meV 
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FIG. 2. Two-triplon contribution to the spectral densities 
of the observable A/"(q) representing the INS on SrCu2(BOa)2 
for two sets of parameters. For clarity, a broadening of 0.02 Ji 
is used. The inaccuracy of the peak positions due to the 
extrapolation is about 2.5%. 



To understand the 2-triplon states throughout the BZ, 
we calculate the spectral density for 150 different mo- 
menta. Color-coding the intensities leads to Fig. 3 where 
we follow the experimentally traced path in dual space. 
The black lines are the most relevant eigen-energies of 
Hi + H2 extracted from the 84x84 matrix representing 
the full matrix of H2 plus a part of Hi at fixed q and 
x. Enlarging this matrix from 84x84 to 112x112 does 
not lead to visible changes. The dashed lines mark the 
lower and upper bound of the 2-particle continuum de- 
rived from the 1-triplon dispersion [9] . The energy range 
depicted is chosen according to a recent high resolution 
INS measurement of SrCu 2 (B0 3 ) 2 [16]. The 1-triplon 
contribution would appear as a sharp, flat and highly in- 
tensive (red) band at about 3 meV from which no new 
insight is gained. 
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Our results compare excellently with the experimental 
data. The main conclusion is that the rather flat bands 
of the 2-triplon states can indeed be understood. Pre- 
viously, experiment [17] and theory [3, 18, 19] found evi- 
dence for significant correlated hopping of two triplons. 
So it came as a surprise that high resolution INS showed 
very flat features only. Previous results were limited in 
resolution (in momentum and in energy [17]), analysed 
only two points of the BZ [3] or were restricted to low 
values of x [18,19]. Fig. 3 shows that there are many 
bound states distributed over a fairly large energy range 
of about 1.5meV. This range corresponds to the previ- 
ous expectation of enhanced correlated hopping. But the 
smoothly connected eigen-energies do not display a sig- 
nificant dependence on momentum. We interpret this 
finding as evidence for level repulsion. Due to the neg- 
ligible 1-triplon kinetic energy there is a relatively large 
number of individual states involved. Their energetic re- 
pulsion renders each individual band very flat. 




q= (2,0) (2.5,0) (2.5,0.5) (2,0) 



K= (0,0) (71/2,71/2) (0,71) (0,0) 

FIG. 3. Color-coded spectral density of M in the en- 
ergy-momentum plane. Intensities as indicated by the scale 
on the right hand side. Parameters are the same as in Fig. 2 

We also computed the energy and momentum inte- 
grated weights. At x = 0.603, we find that about 
50%+25% of the full weight, known from a straightfor- 
ward sum rule, is covered by the 1- and 2-triplon ex- 
citations, respectively. The remaining 25% must be at- 
tributed to higher triplon-excitations. This finding agrees 
nicely with experiment, see e.g. the constant momentum 
scan in Fig. 2(a) of Ref. [17]. 

In conclusion, we like to stress three main results, (i) 
The scenario of strongly dispersing modes in the 2-triplon 
sector cannot be held up. We showed that these modes 
are also rather flat and we argue that this stems from 



level repulsion, (ii) The quantitative agreement with the 
high resolution INS [16] is very good. Not only the over- 
all shape of the structure factor but also prominent de- 
tails like the intensive flat modes at about 5meV and the 
modes just below the continuum at 5.75meV are repro- 
duced. This observation supports our choice for the pa- 
rameters x — 0.603, J\ = 6.16 meV. (hi) Finally, we have 
demonstrated that perturbative CUTs are capable and 
well-suited to quantitatively calculate complex quantities 
like spectral densities also for two-dimensional models. 

We are indebted to N. Aso, K. Kakurai and coworkers 
for making the INS data available to us prior to publica- 
tion. Financial support by the DFG in SFB 608 and in 
SP 1073 is gratefully acknowledged. 
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